Nonaxisymmetric instability of rapidly rotating black hole in five dimensions 
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We present results from numerical solution of Einstein's equation in five dimensions describing 
evolution of rapidly rotating black holes. We show, for the first time, that the rapidly rotating 
black holes in higher dimensions are unstable against nonaxisymmetric deformation; for the five- 
dimensional case, the critical value of spin parameter for onset of the instability is « 0.87. 
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I Introduction: Black holes (BHs) are the most strongly 
self-gravitating objects in nature and also the simplest 
celestial objects, described by a small number of param- 
eters. Therefore, BHs are expected to reflect the natures 
of gravitational theories in the clearest manner. This 
fact motivates extensive studies of BHs not only for four- 
dimensional (4D) spacetime (see, e.g., [l[ for a historical 
review) but also for higher-dimensional one (see, e.g., 
for a review). 

Since a possibility of BH formation in particle ac- 
celerators was pointed out, studies for BHs in higher- 
dimensional spacetimes have been accelerated. If our 
space is a 3-brane in large [|[ or warped Q extra dimen- 
sions, the Planck energy could be of O(TeV) that may be 
accessible with huge particle accelerators like the Large 
Hadron Collider (LHC). In the presence of the extra di- 
mensions, mini BHs may be produced in the accelerators 
and its evidence may be detected. 

A hypothetical phenomenology in the accelerators is 
as follows [!, [f| : In a high-energy particle collision of a 
sufficiently small impact parameter, two particles merge 
to form a distorted rotating BH, and then, it perhaps 
relaxes to a quasistationary state after emission of gravi- 
tational waves. Then, the BH will evaporate by Hawking 
radiation due to a quantum-field-theory effect in a curved 
spacetime. The BH formation and subsequent evolution 
by gravitational radiation can be described by classical 
general relativity [Tl| . Since any approximation breaks 
down for this phase because of its highly nonlinear na- 
ture, numerical simulation in full general relativity is the 
unique approach (see @, H, H, G3 for the 4D case). 

One of the most important issues to be clarified is what 
type of a black object is formed and whether it is stable 
or not. In the 4D case, the formed object has to be a Kerr 
BH because of the uniqueness theorem of BHs (e.g., [l2[ 
for a review), and the numerical analysis of a Kerr BH 
strongly suggests that it should be stable [HI]. These 
facts strongly constrain the possible scenarios for the BH 
formation and subsequent evolution. By contrast, there 
is no uniqueness theorem in higher dimensions: In the 
five-dimensional (5D) case, the black ring solution with 
the ringlike horizon is known in addition to the Myers- 
Perry BH solution with the spherical horizon (but see 
[l4| for uniqueness of 5D black holes with the spherical 
horizon topology). Even if we assume that the black 



rings are unlikely to be formed in two-particle systems, 
the scenario of mini BHs at accelerators is still uncertain 
because there is no proof for the stability of the higher- 
dimensional BHs 3. 

Higher-dimensional BHs with certain parameters are 
known to be unstable against axisymmetric perturba- 
tions. Emparan and Myers suggested that rotat- 
ing BHs with a high spin parameter is unstable for the 
spacetime dimensionality D > 6. The reason is that the 
rapidly rotating BHs have a high degree of ellipticity (i.e., 
the black membrane limit) and such objects are subject 
to the Gregory-Laflamme instability [16j . Very recently, 
Dias et al. indeed showed, by a linear perturbation analy- 
sis, that rapidly rotating BHs for 7 < D < 9 are unstable 
against axisymmetric multiple-ring- like deformation [13] • 

On the other hand, little is known for the stability of 
BHs against the nonaxisymmetric perturbation, and also, 
for the dimensionality D = 5. Emparan and Myers [l5j 
discussed the possibility that the rapidly rotating BHs 
may be unstable also against nonaxisymmetric perturba- 
tion for D > 5 using a thermodynamical argument (i.e., 
by comparing the horizon area of a rotating BH and that 
of two boosted Schwarzschild BHs with the same total 
energy and angular momentum). However, the corre- 
spondence between the thermodynamical and dynamical 
instabilities has not been well established, and thus, a 
rigorous analysis is required in order to clarify whether 
the rapidly rotating BHs are actually unstable against 
nonaxisymmetric perturbation or not. 

In this paper, we tackle this stability issue of rapidly ro- 
tating 5D BHs by fully solving Einstein's equation. The 
merits of this approach are that a wide variety of insta- 
bilities can be clarified with no ambiguity, and that the 
final fate after the onset of the instabilities could be de- 
termined since the amplitude of the perturbation from 
the background BH solution need not be assumed to be 
small. In this letter, we focus on the 5D BHs of only one 
spin parameter because such a BH will be an outcome 
in the particle accelerators. We shall explicitly show, for 
the first time, that rapidly rotating 5D BHs are unstable 
against nonaxisymmetric distortion. 

II Setting and Methodology: We numerically study 
the stability against nonaxisymmetric deformation of a 
5D rotating BH of single spin parameter. Its line element 
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in the Boyer-Lindquist-type coordinates is [18| 

ds 2 = -dt 2 + ^(dt - a sin 2 8dip) 2 + ^-df 2 + Zd6 2 

+ (f 2 + a 2 ) sin 2 6>^ 2 + f 2 cos 2 6»dx 2 , ( 1 ) 

where /i and a are mass and spin parameters, respec- 
tively, and £ := f 2 + a 2 cos 2 6> and A :— f 2 + a 2 — ^. 
r = r/j := (/i— a 2 ) 1 / 2 is the radius of event horizon. In five 
dimensions, q := \a\/ [i 1 / 2 has to be less than the unity 
(i.e., | a. | < /i 1 / 2 ) by the requirement of global hyperbol- 
icity. Note that we adopt the geometric units G = 1 = c 
throughout this letter. 

We evolve this BH adopting the so-called quasi- 
isotropit coordinate r defined by f = r + r 2 /4r. Here, 
the location of the horizon is r — r/j/2. With this trans- 
formation, the t — const, hypersurface becomes spacelike 
everywhere for < r < oo, and the physical singularity 
is not included in the initial surface t = 0. In this co- 
ordinate, the event horizon corresponds to the wormhole 
throat and, the sphere denoted by r = represents space- 
like infinity of another asymptotically flat region and be- 
comes a coordinate singularity. However, the puncture 
approach of numerical relativity (with the appropriate 
choice of the conformal factor and gauge conditions) en- 
ables us to stably follow the evolution of BHs not only 
for 4D spacetimes [22| but also for 5D ones [20| . 

Einstein's equation for a 5D vacuum spacetime are 
solved in the higher-dimensional version of the BSSN 
(Baumgarte-Shapiro-Shibata-Nakamura) formalism [l9l . 
|20| with a fourth-order finite differencing scheme in space 
and time. All the equations are solved in the 4+1 form 
with the Cartesian coordinates (x, y, z,w) where w de- 
notes the coordinate of the extra dimension. The follow- 
ing so-called puncture gauge conditions are adopted for 
the lapse function a and shift vector /?*, 

d t a = -l.haK, (2) 
d t /3 i =r }B B i , d t B l = d t r - iT^B*, (3) 

where B % an auxiliary gauge variable, K the trace part 
of the extrinsic curvature, T l is an auxiliary variable for 
the BSSN formalism, and we choose tjb = 1/3. 

For this BH, rotational motion exists only in the (x, y)- 
plane. Because we consider the nonaxisymmetric stabil- 
ity against a bar deformation, the z- and w-axes direc- 
tions are equivalent. Thus, a rotational Killing vector 
(d/dx)^ is present where \ — tan _1 (Wz). This fact 
motivates us to adopt a cartoon method HH that en- 
ables to solve 4+1 Einstein's evolution equation in the 
3D grid (x,y,z) (see [20j). 

Numerical simulation was performed with two codes. 
One is the code reported in [2(J. This code prepares 
a nonuniform grid for which the region around the BH 
is resolved with a high accuracy and the distant region 
is covered by a sparse grid spacing. The other is the 
SACRA5D code in which an adaptive mesh refinement 
algorithm is implemented. This is the code extended 
from the SACRA code which was originally developed for 



simulations of 4D spacetimes (23|. SACRA5D has been 
tested by solving the problems listed in [2(j. Several 
simulations were performed for both codes and we have 
confirmed that the two codes provide the same conclusion 
concerning the stability of BHs. 

Although the numerical solutions show a behavior of 
convergence with improving grid resolution, the conver- 
gence speed depends strongly on the spin parameter, 
q. In this work, we measured the accuracy by moni- 
toring the area of the apparent horizon by evolving non- 
perturbed rotating BHs and determined the required res- 
olution: We evolved the BHs at least for 50/z 1 / 2 and 
checked required resolution with which the horizon area 
remains approximately constant within 1% error. For 
q < 0.7, we found that the grid spacing near horizon 
Ax w 0.03/i 1 / 2 is small enough for this requirement. 
However, for q > 0.8, the required value for Ax changes 
steeply: For q = 0.8 with Ax/^ 1/2 = 0.015 and 0.0225, 
the errors at t = 50/x 1 / 2 are 0.2% and 1.2%, respec- 
tively. For q = 0.85 and 0.88, Ax/fi 1 ^ 2 should be smaller 
than « 0.01 and 0.008, respectively, to guarantee the 
error within 1% at t = 50fi^ 2 . For q > 0.89, accu- 
rate simulation up to t = 50/i 1 / 2 is not feasible even for 
Ax = 0.008^ 1/2 . Therefore, for this case, the simulations 
were performed up to the time at which numerical error 
of apparent horizon area exceeds 3% error (at t < 50^ 1/2 ) 
choosing Ax = 0.008/i 1 / 2 . Even these short-term simula- 
tions are long enough to show that rapidly rotating BHs 
are unstable. 

To investigate the nonaxisymmetric stability, we ini- 
tially superimpose a small nonaxisymmetric perturba- 
tion on the rotating BH solution. Specifically, we per- 
turb a conformal factor of the 4D space defined by 
W = det(7y)- 1 /4 ( 7ij i s t he 4D space metric) as 

W = W [1 + Afx-^x 2 - y 2 ) exp(-r 2 /2r 2 )], (4) 

where Wq is the nonperturbed solution. In the follow- 
ing, we choose A = 0.005. Simulations were also done 
with values A = 0.01 and 0.02, but the evolution of the 
perturbed part obeys a scaling relation (e.g., W/(WqA) 
behaves approximately in the same manner). Thus, the 
magnitude of A does not change the conclusion in this 
paper, as far as A <C 1. 

During numerical simulation, we monitor two quanti- 
ties for determining the stability against the nonaxisym- 
metric deformation. One is a deformation parameter 
of BH horizon. To define this parameter, we calculate 
the circumferential radii of the apparent horizon along 
several meridians. Specifically, we measure the proper 
length of the meridians (referred to as l v ) for cp = 
(and 7r), 7r/4 (and 5ir/4), tt/2 (and Sir/ 2), and 37r/4 (and 
77r/4), and then, define the deformation parameter 

V= \fik- K/2) 2 + - hir/iY/lo (5) 

The other is the gravitational waveform in a wave zone 
defined along the z-axis by 

ft + =r 3 /V 3/4 (7**-7 w )/2, (6) 
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FIG. 1: Evolution of deformation parameter r\ for a//i 1//2 = 
0.80-0.89. 



where jij is the conformal 4D metric. This quantity is 
regarded as the + mode of gravitational waves. 

In the 4D case, the amplitude of these quantities de- 
creases exponentially with time as far as q < 1, because 
the BHs are stable; the damping rate is primarily deter- 
mined by the fundamental quasinormal modes. However, 
this is not the case in the 5D case (see below). 
Ill Results: Figure 1 plots the evolution of 77 as a func- 
tion of time for q = 0.8-0.89. Irrespective of the values 
of q, 77 quickly increases for t < 10/Lt 1 / 2 . This is simply 
because the perturbation initially given relaxes. Indeed, 
the value of rj after this relaxation phase is of O(10 -3 ), 
which is the same order as the initial perturbation am- 
plitude. However, the evolution after the relaxation de- 
pends strongly on the values of q. For q < 0.86, the value 
of r\ decreases; in particular for q < 0.80, this damp- 
ing occurs very quickly. By contrast, for q > 0.88, the 
perturbation grows exponentially for tj 1 y}l 2 > 25. The 
growth rate does not depend on the initial perturbation 
amplitude. Thus, the BHs of q > 0.88 are dynamically 
unstable against the nonaxisymmetric deformation. For 
q = 0.87, the value of r\ grows slowly with time. This 
suggests that this value is near the critical value for the 
onset of this instability. 

Figure 2 plots h+ and \h+\ as a function of a retarded 
time i rc t =t — r. The amplitude for t > 25/z 1 / 2 exponen- 
tially increases with time for q = 0.89 as in Fig. 1. 

From gravitational waveforms, it is possible to extract 
the frequency of gravitational waves, /. Figure 3 plots a 
characteristic frequency of gravitational waves as a func- 
tion of q. Here, the characteristic frequency is determined 
by performing the Fourier transformation of h + (t) and 
then by identifying its peak. Figure 3 shows that / in- 
creases with q for q > 0.8. This property agrees with that 
of the fundamental quasinormal mode for 4D Kerr BHs 
[li]]. Furthermore, the frequency determined for a = 
gives / w 0.15//^ 1/2 which agrees with that obtained by 
a perturbation analysis [25[. Hence, it is natural to con- 
sider that emitted gravitational waves are associated with 
the quasinormal modes of the BHs. 

A remarkable fact is that the angular velocity defined 
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FIG. 2: h+ and its absolute value as functions of retarded 
time for a/^i 1//2 = 0.85, 0.87, and 0.89 (dashed, long-dashed, 
and solid curves). h+ is extracted for r ~ A where A is a 
gravitational- wave length. 
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FIG. 3: Characteristic frequency of gravitational waves as a 
function of aj p}^ 2 . 



by uj := 7r/ is smaller than that of BHs = aj{r\ + 
a 2 ) = a/fi for a large value of a. Remember that the 
total radiated energy (AE) and angular momentum (A J) 
approximately obey a relation ujAJ — AE, and the first 
law of the BH thermodynamics allows us to obtain the 
variation in the BH area 8 A as kSA/8tt = fl5J — 5E 
where k is the surface gravity of the BH horizon, and 
5E and 5 J are the variation of the energy and angular 
momentum of the BH. If 5E and 6 J are equal to AE and 
A J, respectively, we obtain the relation 



kSA = 87r(f7 - lo)5J. 



(7) 



For fl > u>, S A becomes positive, and thus, the evolution 
by emission of gravitational waves is allowed: A rapidly 
rotating BH in five dimensions may be unstable against 
gravitational radiation reaction (although this is not a 
sufficient condition). This fact seems to be at least a 
part of the reason why rapidly rotating BHs are unstable 
against nonaxisymmetric deformation. 
IV Summary: We showed, for the fist time, that 
rapidly rotating BHs in 5D vacuum spacetime are un- 
stable against nonaxisymmetric deformation. The criti- 
cal value for the spin parameter is « 0.87. The critical 
value we find is close to the value predicted by Emparan 
and Myers, w 0.85 [15j. This suggests that their heuris- 
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tic thermodynamic argument concerning the stability of 
higher-dimensional rotating BHs may be reliable. 

Unfortunately, the present numerical simulation can- 
not clarify the final fate of the unstable BHs, because it 
is not easy to maintain the numerical accuracy to follow 
the unstable BH for a sufficiently long time. To clarify 
the final fate, a simulation with a much better grid reso- 
lution is required. This issue is left for the future study. 

There will be at least two possible fates for an unstable 
BH. One is that the perturbation grows until the unsta- 
ble BH fragments into two BHs [l5|, and the other is 
that the growth of the perturbation saturates at a stage 
when the emission rate of gravitational waves is large 
enough to quickly carry angular momentum of the BHs 
for stabilization. This issue is quite similar to the nonax- 
isymmetric dynamical stability of rotating stars in four 
dimensions: For many cases, rotating stars are dynam- 
ically unstable if the ratio of rotational kinetic energy 
to gravitational potential energy is larger than 0.27 or 
the ratio of the polar axial length to the equatorial axial 
length is smaller than ~ 0.2 [26[. This condition holds 
irrespective of equations of state, as far as the degree of 
differential rotation is not extremely large (e.g., [27] and 
references therein). Dynamically unstable rotating stars 
evolve after the onset of the instability via several mech- 
anisms such as angular momentum transfer and gravi- 
tational radiation reaction. Then, they settle to a new 



stable state. A noteworthy fact is that fragmentation 
rarely occurs for stars (for tori, it may occur). 

For the 5D BHs, the ratio of the meridional circum- 
ferential length to the equatorial one decreases with q, 
and C m /C e = 0.38 for q = 0.87. Thus, a BH is un- 
stable for C m /C e < 0.38. For the 4D case, the mini- 
mum value of this ratio is at most C m /C e = 0.64 (for 
the extreme Kerr BH); all the BHs are not very oblate. 
This may be the reason that BHs are stable in four di- 
mensions. These facts suggest that a sufficiently oblate 
BH with C m /C e < 0.4 is dynamically unstable against 
nonaxisymmetric deformation irrespective of the dimen- 
sionality, as in rapidly rotating stars. Rapidly rotating 
BHs in any higher dimensions can have small values of 
C m /C e [2|, and thus, this issue should be also explored 
by numerical-relativity simulation. 
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